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Abstract: The ATLAS and CMS experiments at the Large Hadron Colhder (LHC) 
may discover the squarks (g) and gluino {g) of the minimal supersymmetric standard 
model (MSSM) in the early stage of the experiments if their masses are lighter 
than 1.5 TeV. In this paper we propose the sub-system variable (771^2'^), which 
is sensitive to the gluino mass when > m^. Using it with the inclusive mT2 
distribution proposed earher, q and g masses can be determined simultaneously in 
the early stage of the experiments. Results of Monte Carlo simulations at sample 
MSSM model points are presented both for signal and background. 
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1. Introduction 

Supersymmetry (SUSY) provides an elegant solution to the hierarchy problem in 
the Standard Model (SM) Higgs sector |T|, ^. It predicts a set of new particles 
containing spin sfermions and spin 1/2 gauginos and higgsinos. If R parity is 
conserved, the lightest supersymmetric particle (LSP), which is often the lightest 
neutralino, is stable and a good dark matter candidate. The thermal relic density of 
the LSP can be consistent with the cold dark matter density of our Universe. 

The ATLAS and CMS experiments at the CERN Large Hadron Collider (LHC) 
may discover the SUSY particles in the early stage of data collection. The missing 
momentum carried by the stable LSP becomes an important signature of the sparticle 
production. Current studies show that an integrated luminosity of around 1 fb~^ is 
enough to find sparticle production if the squark and gluino masses are below 1.5 
TeV and the mass difi^erence between the LSP and squark/gluino is large. 

We do not yet have many clues on the sparticle mass scale, although the current 
measurements of fiavor changing neutral current (FCNC) give stringent constraints 
on the relation among sfermion masses unless they are extremely heavy. Once we have 
seen signs of SUSY at the LHC, we should use direct evidence to determine the SUSY 
particle masses, from which we may determine the sparticle mass relations. Various 
methods have been developed for spaticle mass determination from event kinematics. 
The invariant mass distributions of various exclusive channels are known to be very 
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useful 0, |, H, ^, |, 0. By combining the measured endpoints of the distributions 
of the relatively clean and long cascade decay channels involving charginos (xf)} 
neutrahnos (x^) and sleptons (Z), such as the opposite sign same flavor lepton signal 
arising from q xHl ~^ ~^ XiQ^^^ o^i^ can determine not only the masses of 
the squark and gluino, but also the neutralino and slepton masses arising from their 
cascade decays. The exact relations among momenta of visible particles from a 
cascade decay are also useful |]TD|, |TT|, [T^, ^ . For some cases, the decay kinematics 
can be solved event by event to obtain the sparticle masses in the decay cascade. 

Another important quantity is the 777,^2 variable, which is calculated from two vis- 
ible momenta p^^ig{i = 1, 2) and the missing transverse momentum as in Eq. (|l.l| ) 

BUI, 



= nim 



max {mT(pS, l^?; Mtest), I^?! ^test)}] • (l-l) 



Here, p^^^ii = 1,2) is the sum of momenta of the particles in the visible system i 
which is a set of visible decay products from a parent particle i. M^est is an arbitrary 
mass parameter called the test LSP mass. 

The quantity is bounded above by the mass of the heavier of the initially pro- 
duced sparticles if we set Mtest = "^lsp- This property of the mr2 is useful for 
determining the sparticle masses. For example, for the process pp — >■ gnQR — > jjXiXi 
studied in [0, the events are populated near the 777^2 endpoint, which is very 



clearly visible and coincides with tti^^. 

It was pointed out recently that the endpoint of the 7777-2 distribution 777^2^^ as a 
function of Mtest has a kink at the true LSP mass in the case that the invariant mass 



of the visible system 777vis (which consists of jets and leptons) can range |T^, ^1 



22| , p4| . The kink arises because the derivative of with respect to a test LSP 



mass differs depending on the 7r7vis; while 777!^'^2'^ cannot be above the parent sparticle 
masses max{7rti, 7^72} = M at Mtest = tt^lspi so every trajectory that an event makes 
on the mT2 — Mtest plane passes or goes below the point (7/27-2, Mtest) = {M,mLSp)- 
The LSP mass and gluino mass may be reconstructed from Mtest and the value 
at the kink position. The mass determination has been demonstrated in the four jets 
and channel at a certain MSSM model point in which g decays dominantly via 
g qq —* XiQQ and q is heavy |]21]]. This shows that short hadronic decay chains 
can also contribute to sparticle mass determination. 

To make use of the SUSY events fully in the early stage of the LHC, it is useful 
to define the 7/77-2 variable in an inclusive manner without any specification of decay 
modes. This is because a squark and a gluino may decay into a mode with more 
than two jets in the final state. For example, the decay modes g tixi and tbxf 
have large branching ratios in large mSUGRA parameter regions, because the scalar 
top and the scalar bottom tend to have masses much lighter than the first generation 
squark masses. 
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In the previous paper, we therefore define an inclusive mT2 variable using a 



hemisphere method [^. The inclusive 112x2 is defined in two steps. In the first step, 
we divide jets in each event into two hemispheres ||25|, This is normally done by 
associating the jets with two leading axes which are initially taken as the two leading 
jets in the event. The sum of the jet momenta assigned in a hemisphere is called 
a hemisphere momentum p[,emi(^ = 1)2). In the next step, a stransverse mass 
is calculated using Eq.(|l.l|) with p^ij^i taken as p^^. The inclusive 171x2 as defined 
above carries the information on the parent sparticle masses max{mi,m2}, if the 
hemisphere algorithm groups the decay products from the particle 1 and 2 into two 
different hemispheres correctly. It is shown that a parent squark mass rriq can be 
obtained from the case of rriq > nig. Moreover, the as a function of a 

test LSP mass still has a kink at the true LSP mass. The inclusive 171x2 distribution 



is also useful for discriminating model parameters and discussed extensively in [27 



In this paper, we propose a "sub-system mx2" , Tn^2 ■ defined as an inclusive 
mx2 variable, but the highest px jet is removed before the hemisphere reconstruction. 
The definition is inspired by an observation that the squark decays via q ^ g ot 
q produce a high px jet if rriq is sufficiently larger than rrig, m^± and 

m^o. In the case that the jet from the squark decay is identified, the remaining 
system is either gluino-gluino or gluino-neutralino/chargino for qg co-production, so 
m^2^(end) = rrig. By studying several model points we provide convincing cases that 
both rriq and rrig are estimated using ni^2 and m^g^'*''^*^. We also calculate background 
distributions coming from the productions tt-t-ra jets, W + n iets and Z + n jets using 
ALPGEN [g8|, m with MLM matching. We find that the signal to noise ratio {S/N) 
is large especially for the events near the mx2 endpoint, which are most sensitive to 
squark and gluino masses. 

The importance of matrix element (ME) corrections to SUSY processes have 



been emphasized recently [pQ, 31, 32 1. We provide an estimate of the size of the 



matrix element correction to the signal mx2 distributions using MadGraph |^3|]. We 



find that the signal mx2 distributions are not significantly modified by the SUSY 
matrix element corrections near the endpoint of the mx2 distributions. 

This paper is organized as follows. In Section |^, we describe the sub-system 
mx2- We show parton level mx2 and m^2 distributions, and discuss reconstruction 
efficiencies of the SUSY decay cascades using a hemisphere algorithm at our sample 
model point. The jet level distributions using HERWIG ||3^ with simple detector simu- 



lator AcerDET are given in Section ^. We study the SM background distributions 
generated by ALPGEN in Section |^. Section | is devoted to the conclusions. 



In Appendix |A.1| , we describe the condition that m^2 coincides with the mass 



of the heavier of the initially produced squarticles. The effects of the matrix element 
corrections to mx2 distributions are studied in Appendix [A. 2. 



3 



2. The sub-system inclusive mT2 (^r2^) - Parton level analysis 

At the LHC, squarks and gluinos are copiously produced via qq, qg and gg production 
processes. Each of them decays into visible objects and a LSP. If the visible systems 
are correctly grouped, the inclusive 771^2 with the correct defined as Eq. (|l.l| ) can 
be calculated. In that case, the important property is^ 

m^2 = max{mi, 777,2}. (2-1) 

Here, denotes the endpoint of the mT2 distribution and mi and m2 denote the 
masses of the produced parent particles. In this section, the test mass is taken as 
the true LSP mass (Mtest = "^lsp)- 

In the case of rriq ^ nig, squark-gluino and gluino-gluino production events 
are dominant SUSY production processes. They give the different mT2 endpoints: 
rriq and rrig. For the squark-gluino production events, a squark decays dominantly 
into a gluino (or another lighter sparticle) and a jet. If we can identify the jet, 
all other elements of the system make up a sub-system that may be considered as 
a gluino-gluino (or gluino-the other sparticle) system. We introduce the variable 
171^2 (sub-system mT2), which is mx2 calculated for the sub-system. The missing 
transverse momentum is taken as the same as for the whole system since the sum of 
the two LSP momenta is required for the calculation of mT2- The expected endpoint 
of m|^2^ is rrig. 

Practically, we define the sub-system as the system with the highest pr jet re- 
moved. If the highest pt jet is from a decay chain of a squark, the endpoint of m?^2 
is expected as, 

777,^2*''^'^'^ = min{mi, m2}. (2.2) 

We now show parton level mT2 and m?^2 distributions at our sample model 
points. Here, we take the model points "a - f" listed in Table |1| with the GUT scale 
gaugino mass M1/2 = 300 GeV and tan/3 = 10. The GUT scale sfermion mass mo 
is 600 GeV at point f and 1400 GeV at point a. The gluino masses at these points 
are approximately the same, while squark masses range from ~ 900 GeV (at point 
f) to ~ 1500 GeV (at point a). The GUT scale Higgs masses are tuned so that the 
/X parameters are small ~ 180 GeV. The relation fi ~ Mi results in a thermal relic 
density of the LSP that is consistent with the observed cold dark matter density of 
our universe 



36| , [37| , p8[ . Some of the branching ratios of the 1st generation squarks 



are given in Table 0. We calculate the masses and the branching ratios at these 



^The endpoint for qg production events is given by inax.{mq,mg} unless the LSP mass is too 
close to rrig, which is satisfied in the typical mass spectrum. More details are discussed in Appendix 



A.l 



^The choice of the fi parameter does not affect mT2 distributions discussed in this paper. 
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ITLq 


u 


rriq 


9 


-LjO-T 




a 


1400 


-1400 


1516 


795.7 


107.9 


180 


b 


1200 


-1200 


1342 


785.0 


107.4 


180 


c 


1100 


-1100 


1257 


779.5 


107.1 


180 


d 


1000 


-1000 


1175 


773.2 


106.8 


180 


e 


820 


-750 


1035 


761.7 


106.1 


180 


f 


600 


-650 


881.0 


745.4 


107.8 


190 



Table 1: Some of the mass parameters of our model points. We take the scalar masses of 
sfermions and gaugino masses to be universal. We tune the higgsino mass parameter // by 
allowing non-universal GUT scale Higgs masses parameters so that Qh'^ ~ 0.1. All mass 
parameters are given in GeV. 



point 




Br{uL — > gu) 


Br{uR gu) 


(T(SUSY)(pb) 


a{q) (pb) 


a 


1516 


0.71 


0.93 


4.91 


0.46 


b 


1342 


0.68 


0.92 


5.35 


0.79 


c 


1257 


0.66 


0.91 


5.84 


1.07 


d 


1175 


0.62 


0.90 


6.15 


1.40 


e 


1035 


0.53 


0.96 


7.31 


2.36 


f 


881 


0.31 


0.71 


9.49 


4.34 



Table 2: Some relevant branching ratios of squarks are calculated using ISAJET. The total 
SUSY production cross section and the cross section involving at least one first generation 
squark estimated using HERWIG are also given. 



model points using ISAJET and the mass parameters are interfaced to HERWIG 
using ISAWIG |^0|. The cross sections are calculated using HERWIG. Note that the 
values of mo and M1/2 in Table ^are within the discovery reach in mSUGRA for the 
ATLAS and CMS experiments at / dtL = Ifh'K 

The squark production cross section reduces quickly with increasing first gener- 
ation squark masses. The total SUSY production cross section varies more than a 
factor of two from point a to point f. The difference comes mostly from the decrease 
of CT{qg). In particular, the production cross section involving at least one first gen- 
eration squark is only 0.46 pb at point a and 4.34 pb at point f. The gluino-gluino 
production cross section also becomes reduced because t-channel squark exchange is 
suppressed. Chargino and neutralino production is important at point a. 

The squark decays dominantly into gluino and a jet (See Table 2). The squark 
branching ratio into the gluino is dominant except at point f. For points a to c, 
the mass difference between squark and gluino is significantly larger than half of the 
gluino mass. Therefore, a jet from a squark decay is likely to be the highest pt jet 
in the event. Jets from the other squark decay modes such as q ^ Xij and x^j have 
Pt which is even higher than that of g ^ gj on average. 
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To define the inclusive mT2 and m?p2 distributions, we group jets in an event into 
two "visible objects". For this purpose, we adopt the hemisphere method in Refs 



25, Eq]. For each event, two hemispheres are defined and high jets are assigned 



to one of the hemispheres as follows: 

1. Each hemisphere is defined by an axis Phemi ~ 1) 2), which is the sum of the 
momenta of the selected high objects belonging to the hemisphere i. 

2. A high pj- object k belonging to a hemisphere i satisfies the following conditions: 

d{Pk,PhLi) < diPk,PhLi), (2-3) 



where the function d is defined by 



diPk.P^U = (4L - IpILI cosg.,) f^-' (2.4) 

l^hemi + ^k) 

Here dik is the angle between p^l^^ and p^. 

For our 171x2 analysis in this paper, the selected objects are jets withp^-j > 50 GeV 
and Ir^jl < 3. We do not include the jets with px < 50 GeV nor |?7j| > 3 to avoid 
contaminations from initial state radiations. For m?^2> we do not include the highest 
Pt jet in the selected objects. 



To find the Phemi' adopted the algorithm discussed in Refs. p5|l : 



1. We first take the highest px object with momentum pi and the object i which 
has the largest pxi x Ai?(l,z), where AR{i,j) = 

2. We regard pi and Pi as two seeds of the initial hemisphere axes, and assign the 
other objects to one of the two axes. 

3. We recalculate the hemisphere axes. We perform iterations until the assignment 
converges. Once Phemi's are determined, one can calculate mT2 by using Eq.|LT| 
with taking phemi as p^is- 



In this section, we study parton level events. The momenta of quarks and gluons 
from sparticle decays are extracted from HERWIG event records, and only q-q^*\ (f*^-g 
and g-g productions are included in the figures. When a sparticle decays into W^,Z^, 
and t, we further follow their decays. Note that each parton is in general off-shell 
when they are created from a sparticle decay, and we do not follow parton shower 
evolutions after that. We do not include partons from initial state radiations. These 
effects will be taken into account in particle level MC simulations in the next section. 

For the calculation of the m™2^, we remove the highest pt jet before the hemi- 
sphere assignment to obtain the Phemi's. As an alternative definition, we can remove 
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500 1000 1500 [GeV] 500 1000 1500 [GeV] 



Figure 1: Various parton level m|V|' distributions at point b. The left plot, a solid 
line: the m?^2 distribution. The right plot, a solid line: a distribution, which is a 
mj'2 variables but calculated after subtracting the highest pT object from the hemisphere 
momenta. Dotted lines, m™2'^(true) distributions, which uses generator information for the 
hemisphere assignment. Arbitrary normalizations are used for the y-axes. See text for 
details. 



the highest px jet from the Phemi's after the hemisphere assignment, and denotes 
this alternative in the following. 

We now compare ^nd m^2^ at point b in Fig. |l[ In the left plot, we show 
the m?^2 distribution in the solid line. In the right plot, the solid line shows the mT2 
distribution. In each plot, the dotted line shows the 'true distribution' m|if'(true), 
in which the p^^ consists of the momenta of decay products from a parent particle 
i except for the highest p^ parton using the generator information. This is an ideal 
distribution when the assignment of the visible systems is perfect. Note that the 
highest pt jet is not always from a q decay. Even in the distribution of m^2''(true), 
two endpoints can be seen, the lower is at the gluino mass and the higher is at the 
squark mass. 

The endpoint at the gluino mass is more clearly visible for the m?p^ than for 
the rh^2 distribution. The improvement in m^2 distribution may be explained as 
follows. At point b, a parton from q ^ qg has a large open angle to the gluino decay 
products on average. The event effectively has three axes: the two momenta of the 
two gluino decay products and the momentum of the extra parton from squark decay. 
The assumption of the hemisphere algorithm that events must have two axes may 
lead to an incorrect hemisphere assignment. Removing the highest px parton before 
the hemisphere assignment therefore makes the kinetic endpoints more visible. 

The successful endpoint reconstruction shows that the hemisphere algorithm 
reconstructs a total visible momentum of a squark/gluino decay more or less correctly. 
One can check this explicitly by counting the number of partons assigned to an 
incorrect hemisphere. (Fig. |]). The solid (dashed) histograms correspond to the 
distributions of the number of mis-reconstructed partons for the case that the highest 
Pt parton is removed before (after) the hemisphere assignment. The improvement 



-7- 



woo 

800 
600 
400 
200 




Point a 



0123456789 



1600 
1400 
1200 
1000 
800 
600 
400 
200 




Point b 



0123456789 



3000 




2500 




2000 




1500 




1000 




500 







Ml 






Point e 



0123456789 



Figure 2: The distributions of the number of partons assigned to an incorrect hemisphere 
at point a (left), b (middle) and e (right). In all figures, solid (dashed) histograms cor- 
respond to the case that the highest px parton is removed before (after) the hemisphere 
assignment. Arbitrary normalizations are used for the y-axes. See text for details. 




1500[GeV] 



500 1000 1500[GeV] 



Figure 3: The m?p2 distributions (solid) and the m™^(true) distributions (dashed) for 
squark-gluino co-production events at points d and e. Arbitrary normalizations are used 
for the y-axes. 



achieved by removing the highest pr parton before the hemisphere assignment is 
clearly seen at point a (the left plot). At this point, nig = 1516 GeV and = 
796 GeV, so the parton from the squark decay should have pt of the order of several 
hundred GeV. We also see mild improvement at point b (the middle plot). At point 
e (the right plot), the squark and gluino masses are close, (m^ — mg)/mg = 0.36. In 
this case, removing the highest pt jet before the hemisphere assignment leads to the 
slightly worse reconstruction efficiency. The number of mis-reconstructed partons is 
either or one for more than half of the events in Fig. ^. 

If the highest px parton does not arise from q decay, the gluino endpoint cannot 
be reconstructed even for squark-gluino production events. The probability strongly 
depends on the model parameters. In Fig. ^ we show the distributions for 

only squark-gluino co-production events at points d (left) and e (right). The gluino 
endpoint can be seen around 750 GeV from the m?^2 distribution at point d, which 
is close to that of the m|!if'(true) distribution shown in the dotted line. However, at 
point e, even in the m^2^(true) distribution we cannot see the clear structure at the 
gluino mass. 
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Figure 4: The parton level mT2 distributions at points a, b, d and e. The solid lines 
correspond to the distributions using the hemisphere algorithm, while the dotted lines 
correspond to those with correct parton assignments obtained by using the generator in- 
formation. Arbitrary normalizations are used for the y-axes. 



The difference between the m?^2 ^^"^ m™2''(true) distributions at points d and e 
may be explained as follows. At point d, (m^ — mg)/m^ = 0.52, and the energy of the 
parton from the squark decay is bigger than that from the gluino on average. This is 
why m™2^ shows clear gluino endpoints at point d. In contrast, {rriq — rrig) / nig = 0.36 
at point e. The rriq — rrtg is not large enough, and it is not likely that the parton from 
q —>■ gq has significantly high compared with those coming from the gluino decays. 
This is why the m^2^(true) distribution does not show the endpoint at the gluino 
mass, it could be a problem to extract the gluino mass from the m?p^ distribution. 
However, the m?^2 distribution of qg production is significantly smeared towards the 
lower mT2 value. In the actual situation, the contribution from gluino-gluino pair 
productions would be added, and the distribution would have the endpoint 

at the gluino mass. We will see in the next section that the contamination from 
squark-gluino production is not serious. 

For completeness, we show the parton level mT2 distributions at points a, b, d 
and e in Fig. ^ to emphasize the difference between and distributions. The 
solid histograms are the mT2 distributions using the hemisphere algorithm, while the 
dotted histograms correspond to the mT2(true) distributions which are obtained by 
assigning the partons arising from a parent particle i to hemisphere i using generator 
information. At points a, b and d, the mT2(true) distribution has two peaks. The 
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Figure 5: The inclusive m^a destributions at points a b, d and f (from left to right). The 
up squark mass is 1516, 1342, 1175, 881 GeV from left to right, respectively. Unit of y-axes 
is events/bin/1 fb"-*^. 

peak at the lower mT2 value comes from gluino-gluino production, while the peak at 
higher mT2 corresponds to the squark-gluino and squark-squark productions. The 
endpoint of the distributions coincides with squark mass. The double peak structure 
cannot be seen in the distributions of mT2i but the endpoints are the same as that 
of mT2(true). 

The slope of the distribution near the endpoint becomes flatter with increasing 
squark mass as can be seen from the distributions at points a, b, and d. In particular, 
the existence of a high pt parton from squark decay leads to some confusions in the 
hemisphere algorithm at point a, and a careful study of the distribution would be 
required to extract the squark mass from the fit. The peak of the distribution 
coincides with the lower 771^2 (true) peak. The events near the peak come from gluino 
pair productions at points a, b, and d. In principle, the position of the peak contains 
the gluino mass information. However, this is not easy to observe because the SM 
background may also be large in this region. At point e, although the endpoints of 
the mT2 and mT2(true) distributions are consistent, the squark and gluino masses 
are too close for the two peak structure to be seen. 



3. The MC simulation of the signal 

We have shown that the endpoints of 771^2 and distributions carry the infor- 
mation on squark and gluino masses using parton level events. In this section we 

study the events produced by a parton shower Monte Carlo HERWIG (in the particle 
level) with a detector simulator AcerDET under the set of cuts to reduce the standard 
model backgrounds. The simple snowmass cone algorithm implemented in AcerDET 
is used for finding jets and we set the cone size R = 0.4. 
We apply the following cuts to the events. 

• Jet Pt cuts: ^50 = A'jet(Pr > 50 GeV) > 4, riioo = N^etiPr > 100 GeV) > 1. 
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Figure 6: The fitted mT2 and endpoints (solid lines) and rriq and nig (dashed lines) 
at each model point. The bars show the size of statistical errors for 50,000 SUSY events. 

• Meff(= Ej^^^>50GevPT + $t) > 500 GeV 

• Transverse sphericity: St > 0.2. 

• Missing Transverse momentum: > 200 GeV, J^t > 0.2Mefj. 

• No isolated lepton with px > 20 GeV. 

These cuts are similar to the standard SUSY cuts in the ATLAS TDR 0, except 
for our tighter cut. We veto events with isolated leptons because a hard lepton 
might be assoicated with a hard neutrino. If there is a hard neutrino in an event, 
PTmiss of the event may not be the sum of the transverse momenta of LSPs. In that 
case, the endpoint of the 771x2 distribution might be smeared. 

We first show the mT2 distributions for our model points. The mT2 distributions 
for mtest = 10 GeV under the SUSY cuts are shown in Fig. For each point, we 
have generated 50,000 SUSY events and the distribution is scaled to correspond to 
J dtC = 1 fb~^ of luminosity. 

The endpoint of the mT2 distribution is roughly at ~ nig. We fit the distributions 
to linear functions 

/(m) = a(m — + c (for m < 171^2) 

= b{m - m^^^) + c (for m > m^^^), (3.1) 

and the fitted m^2 values are shown in Fig. p. Here, the statistical errors shown 
in bars correspond to 50,000 total SUSY events. The obtained m^2 rrig are 

^We set mtost small as we do not know the LSP mass initially. 
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Figure 7: The pT distributions of the highest pT jet (left) among the jets with |r/| < 3, and 
the distributions of the m™|' (right) at point a. The dashed hues show the contributions 
from the events with naoo = 1 and the dotted hues show the contributions from the events 
with nq = l. The standard SUSY cuts are apphed for all plots, and mtest = 10 GeV. Unit 
of y-axes is events/bin/lfb~^. 



consistent except at points a and f. For point f, the squark and gluino masses are 
too close, and it is natural that the endpoint fall at weighted mean of gluino and 
squark masses. For point a, due to the very large mass difference between squark and 
gluino, the hemisphere method involving the highest pt jet does not work perfectly. 

Note that there is some ambiguity in choosing a fitting region. For example, 
for point a, the distribution consists of the two components, one arising from the 
gluino-gluino production with the endpoint around 800 GcV and the other from the 
squark-gluino production with the endpoint around 1400 GcV. We fit the distribution 
above mT2 > 1000 GeV for this point. If we did the same fit at point f (the right 
plot), we might fit the mis-reconstructed tail of the events and therefore might obtain 
the endpoint at 1150 GeV. This suggests that the region of the fit must be chosen 
carefully. In particular, the events near the fitted endpoint must make up a sizable 
fraction of the total events. For points b, d and f, we first fit the region from the 
mT2 slightly above the peak position of the distribution up to the highest bin with 
enough statistics (> 10 events/bin). We then increase the lower limit until we obtain 
a small Ax^. The Ax^/n.d.f is less than 1 except at points c and e, and all fits 
satisfy Ax^/n.d.f < 2. 

We now demonstrate the gluino mass determination using the endpoint of the 
mj^2 distribution. Here we must pay some attention to reduce the contributions from 
the squark-squark pair productions, which give the endpoints of the distribution 
as rriq. This contribution smears the endpoint at the gluino mass. It is important to 
find the cuts to reduce the events. 

We find that the cut on the number of high px jets above a certain threshold is 
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Figure 9: The destributions at points a to f normalized for 1 fb~ . The input gluino 
masses differ shghtly among the model points as we fix gaugino mass at the GUT scale. 
The top figures correspond to points a, b and c from left to right, and rrig = 796,785, 
and 780 GeV respectively. The bottom figures correspond to points d, e and f, and nig = 
773, 762 and 745 GeV. The squark mass is 1516 GeV at point a and 881 GeV at point f. 



useful to reduce the contamination, becuase the squark decay tends to give high pt 
jets, as we have discussed earUer. To see this, we first show the distributions of the 
highest Pt jet at point a for 50,000 generated events in Fig. |^ The sohd hues show 
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the pT{ist) distribution, where p{lst) is the momentum of the highest pt jet among 
the jets with \r]\ < 3. The dashed hues show the contribution from the events with 
77-300 = NjetiPT > 300 GeV, Ir^l < 3) = 1 and the dotted hues show the contributions 
of the events with rig = 1, where Ug is the number of primary produced 1st generation 
squarks of the events. The standard SUSY cuts are apphed to the events. We can 
see that most of the events with pxilst) > 300 GeV satisfy ^300 = 1 and they mostly 
come from squark-gluino productions. Therefore, if 71300 ^ 2, they are hkely come 
from squark-squark pair production events. 

Based on the above observation, we calculate the distribution only for the 
events which have only one or zero high px jet above a certain p-r threshold. The 
actual value of the pr cut should be chosen based on the signal distribution. For 
our model points, we take the cut 7^300 < 1- We do not include the events with 
%oo ^ 2, because our MC simulations show that they mostly come from the squark 
pair production. In the right figure, we show 777^2'' distributions for the events 71300 = 1 
at point a. The dotted line shows the distribution with 71300 = 1 and rig = 1. The 
dashed line is the distribution with 111x2 > 800 GeV, 71300 = 1 and rig = 1. All 
distributions show the m?^2 endpoint close to the gluino mass value ~ 800 GeV, 
which is expected from the parton level analysis. 

Fig. ^ shows the same distributions at point f. Events from squark-gluino co- 
production still dominate the events with px > 300 GeV, and a significant fraction 
of the events satisfy 71300 = 1- The events near the m?^2 endpoint mostly come from 
squark-gluino production. The endpoint of the distribution ~ 750 GeV is consistent 
with the gluino mass. 

The distributions at points a to f for 1 fb~^ of integrated luminosity are 
shown in Fig. Here we require 71300 < 1; therefore, the distributions now include 
significant events from gluino-gluino production unlike the previous plots. We have 
seen that the mT2 distribution changes significantly among points a to f. The m|^2^ 
distributions are, by contrast, similar. This is because the 771^2^ endpoints must be 
very close to the true gluino mass nig ~ 750 GeV (up to the difference of the test 
LSP mass from the true LSP mass). This is also seen in Fig. ^, where the value of 
the fitted m?^2 endpoint is shown together with the gluino mass for each point. 

4. Background 111x2 and mf^2^ distributions 

The Standard Model background to the SUSY processes has been studied by ATLAS 
and CMS groups extensively. The ratio /Mes gives a good discrimination between 
the SUSY signal and the background. In the previous section we required I^t/Mss > 
0.2 in addition to M^s > 500 GeV and > 200 GeV. 

The production cross section of the SM background is huge compared with the 
typical signal cross section. To measure the endpoint of the signal 771^2 and m?^2 
distributions, the signal to noise ratio (S/N) must be sufficiently small near the 
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ti W Z total 


mT2 > 500 GeV 


77.7 104.9 107.0 289.6 
38.4 44.8 39.9 123.1 


mT2 > 700 GeV 


20.3 24.4 23.4 68.2 
10.0 12.0 10.4 32.4 


mf^ > 300 GeV 
mf^ > 500 GeV 


90.3 80.4 82.2 252.9 
44.2 38.5 31.7 113.1 
11.1 6.9 6.1 24.0 
8.1 4.9 3.8 16.8 


luminosity 


13.1 fb^i 13.5 fb-i 19.1 fb-i 



Table 3: Number of SM background events per Ifb^^. For each row, upper (lower) 
numbers correspond to the events without (with) a cut on the hemisphere masses, m-hemi > 
200 GeV. The last row shows the number of generated events for this study in terms of the 
corresponding integrated luminosity. 



endpoint. The SM backgrounds in the 0-lepton channel after the standard SUSY 
cuts come from the four different sources: ti, W^, Z° productions with multiple jets, 
and QCD multi-jet processes. Bottom quark productions and the mis-measurements 
of particle energies can give the missing energy to QCD multi-jet processes. It is 
difficult to estimate the QCD background without knowing detector performances 
in detail. We therefore do not attempt to do so in this paper. In recent ATLAS 



and CMS studies |T^, the four channels contribute to the background at roughly 
the same order of magnitude after the cuts to reduce the SM backgrounds, although 
QCD background decreases much faster with increasing Meg. 

The source of missing Et for the processes tt, Z^, and + n jets is primarily 
escaping neutrinos, and missing Et arising from energy mis-measurements is less 
important. We generate these events using ALPGEN P5| , and parton shower and 
initial state radiations are estimated by interfacing the parton level events to HERWIG. 
We generate Z^[-^ uv) + n jets for n < 5 , W^i^-^ lu) + n jets (^ < 4) , and ti + n 
jets {n < 2) , so that tree level lepton events have at least 4 or 5 jets including r 
jets. We require minimum parton separation ARjj > 0.6^, and place a cut on the 
forward parton of |?7| < 5. The events are then matched so that there is no double 
counting between parton shower and hard partons by using the MLM matching 
scheme provided by ALPGEN. In this scheme, we generate the processes with up to 
nmax parton. The events from the processes with n partons (n < nmax) are accepted 
only if jets and partons match (njet = n), while the events from the processes with 



^The jet cone size for AcerDET jet reconstruction is set to i? = 0.4. This means that centers of 
two well separated jets has AR > 0.8. We therefore require ARjj to be slightly lower than that. 
This is sufficient for our purpose as we are working on inclusive signatures. Reducing the ARjj cut 
to less than 0.6 results in unnecessary inefficiency to the event generation. 
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'^max partons are accepted if n^et > n. In order to reduce the number of produced 
events while keeping enough statistics for the kinematical region we are interested 
in, we require Xlparton > 400 GeV for W+ n jets, I^t > 150 GeV for Z + n jets, 
l^parton-^r > 500 GcV for ti+ n jets^. The effect of additional jets on the signal 
distributions is small and discussed in Appendix [A.2 



AcerDET performs Gaussian smearing for jet momenta, the missing momentum, 
and isolated lepton momenta. It does not contain various potentially important 
instrumental effects, such as non-Gaussian tails of the energy smearing and lepton 
inefficiencies. Therefore, our background estimate is given in this paper for illustra- 
tive purpose, and more realistic estimates must be performed by the experimental 
groups. 

Keeping this in mind. Table ^ summarizes results of our event generations. The 
number of background events for J dtC = 1 fb~^ under various cuts are given. The 
bottom row shows the corresponding luminosities we have generated for the back- 
ground processes. We apply the SUSY cuts given in Section 4. In addition, we 
require u^qq < 1 for m^2^ distribution. We do not include K factors, as the corre- 
sponding higher order QCD corrections are not available. Note that K factors of 
tt production and SUSY production tend to cancel partially. For each row, upper 
(lower) numbers correspond to the events without (with) a cut on the hemisphere 
masses, mhemi > 200 GeV. The background with the hemisphere mass cut is reduced 
by more than a factor of 2. This suggests that the background events are dominated 
by the configurations that a few jets are either soft or colinear to leading hard jets 
and therefore the masses of the hemispheres are small. The background distributions 
will be studied in detail elsewhere. 

Fig. [l^ shows the distribution of background, together with the signal distribu- 
tion at points f (the top figures) and d (the bottom figures). These distributions 
are without hemisphere mass cuts. The signal is larger than the background above 
mT2 > 600 (750) GeV at points f (d) for the mT2 distribution, which is much smaller 
than expected = = GeV ( 1175 GeV). The endpoints of the signal mT2 
and distributions may be extracted as a kink in the total distribution in this 
case. The signal and background distributions of ^ilso shown in the right 

plots. Again, the level of the background is small near the endpoint. 

We also show the same distribution at point b in Fig. [Tl|. The signal cross section 
involving q production is reduced by a factor of 1/5 from that at point f (See Table 
H) . The S/ N above mT2 > 700 GeV is now ~ 1 and the cross point of the signal and 
the background is at mT2 = 800 GeV. By applying the hemisphere mass cut, we can 
reduce the background significantly. The improvement of S/ N near the endpoint can 
be seen by comparing the top and bottom figures without /with the hemisphere mass 



^The conditions of the generations for the different processes are not the same. However, these 
conditions are loose enough so that there is no effect of the generation cuts after our standard SUSY 
cuts. 
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Figure 10: Signal and background distributions at point f (the top figures) and at point 
d (the bottom figures). Unit of y-axes is events/bin/1 fb~^. 



cut. It is important to reduce the background to measure the squark and gluino 
masses near the discovery regions. 



5. Conclusions 

The ATLAS and CMS experiments at the LHC can discover squarks and gluinos in 
the MSSM with masses less than 1.5 TeV at the early stage of the experiment with 
luminosity around J dtC = lfb~^. Developing a reliable method of estimating squark 
and gluino masses with the discovery is an important step to study supersymmetry 
at the LHC. 

For this purpose we cannot rely on the clean golden channels such as 
jets, becuase they tend to have small branching ratios and are sensitive to the model 



parameters. In a previous paper [^], we defined an inclusive mT2 variable. This 
variable can be calculated for any event with jets and missing transverse energy. 
It is calculated in two steps; we first define the two hemisphere axes by assigning 
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Figure 11: The signal (solid) and background (dashed) distributions at point b. The 

1(2) 

bottom figures show the distributions with the hemisphere mass cut m^^^^ > 200 GeV. 
Unit of y-SiX.es is events/bin/1 fb"^. 



particles into the two leading jets of the events, then, the mT2 variable is caluculated 
from the two hemisphere momenta and missing transverse energy. We pointed out 
that the endpoint of the mr2 distribution is sensitive to the squark mass for the case 
rriq > rrig. 

In this paper, we define a "sub-system" mr2, m?p2- This is an mT2 variable 
calculated without including the highest pr jets for the hemisphere assignments and 
772^2 calculation. In the case that rriq > rrig and the other sparticles are lighter, the 
endpoint of distribution gives us information on rriq. In this paper, we show 
convincing evidence for sample model points within the reach for J dtL = 1 fb^^. 

We also provide various parton level checks on the hemisphere algorithm. We 
estimate background distributions arising from ti+n jets, n jets and + 
n jets using ALPGEN and find out that S/N ratio is large for the events near the 
mT2 endpoints at our sample points. In the Appendix, we also provide a study of 
SUSY+n jet distributions using MadGraph/MadEvent, and find that the endpoint is 
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stable with the ME corrections. 
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A. Appendix 

A.l The endpoint for squark-gluino production events 

In this Appendix, we show the condition for which the endpoint of the ideal 
distribution for the squark-gluino production events coincides with the squark mass 

at Mtest = "^LSP- 

The squark-gluino mT2 is calculated by minimizing max{m® , m^'' } under the 
condition that the sum of transverse test momenta of two LSP is equal to the I^t- 
It is known that the transverse mass {i = q, g) as a function of the test LSP 
momentum has the global minimum, which is called the unconstrained minimum 
(UCM) []T3| . There are cases where mT2 is given by the unconstrained minimum of 
the transverse mass on one side (m^^)ucM- This situation occurs when m^'' on the 
other side for the test LSP momentum which gives the (^7i^^)ucM is smaller than 
(m^ )ucM- 



In Ref.[18|, it is shown that the UCM of the squark system )ucm) is given 

by 

rf)uCM=mi^+ Mtest, (A.l) 

where m^f^ is the invariant mass of the visible particles from the squark decay. The 
maximum of (m^^)ucM is, therefore, given by substituting the maximum of m^f^ into 
Eq. |A.l| . The maximum of the is given by 



(m!|))-- = m,--mLSP, (A.2) 

if the LSP from squark decay can be at rest in the squark rest frame. In this 
case, the maximum of the UCM of the squark system can reach the squark mass at 
Mtest = ''^LSP, and the maximum of the squark-gluino mT2 is identical to the squark 
mass at Mtest = "^lsp- 

We now consider the condition that the LSP can be at rest in the squark rest 
frame. In the following discussion, we concentrate on the case that the squark decays 
into the gluino and a jet. The gluino from the squark subsequently decays into the 
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visible objects and the LSP (See Fig. 12). The LSP momentum in the gluino rest 
frame (plsp) depends on the invariant mass of the visible objects (m^f]) as 



IPlspI = ^ V^I - 2^i(("^3])^ + Ksp) + iim^r - mU^r ■ (A.3) 

If the LSP is produced in the opposite direction from the gluino momentum and the 
gluino velocity is not too large, the LSP can be at rest in the squark rest frame for 
a suitable value of the invariant mass of the visible objects (m^Jil)- In this situation, 
the LSP momentum in the squark rest frame (plsp) is obtained by the Lorentz boost 
of the PLgp as 

Plsp = 7gH/3g\Ksp + IPlspI) , (A.4) 

where -^lsp is the energy of the LSP in the gluino rest frame, and the Lorentz boost 
factors Pg and jg are given by 

By solving the equation plsp = 0, we obtain 

(m^j) =m^(l-— — )(1 • (A.6) 

iilq g 

Note that if m^mLSP > "n^^, the equa- p. 
tion plsp = does not have any so- ("^i^^'" '* ^ ^ ^ 



VIS/ 



lution for positive w^^fg. In this case, visible -* *- LSP 

the LSP cannot be at rest in the squark ^.^^^^ ^ Kinematical configuration of 
rest frame, and (m^^)ucM is less than gquark decay 
the squark mass. Even if the equation 

Plsp = has a solution for positive m^f], there are cases where m^f] has a non- 
vanishing kinematical lower bound due to a heavy standard model particle, such as 
t, Z and W. If the lower bound is smaller than the solution m^f], (m®)ucM cannot 
reach the squark mass. For our model points, the solution (|A.6|) is m^f^/rrig = 0.83 
for point a, and fh^^l/mg = 0.85 for point f. On the other hand, the kinematically 
allowed range of the visible invariant mass is roug hly 0.50 < m^£/mg < 0.86 for 
point a, and 0.55 < m^^/rrig < 0.86 for point f. Therefore, the endpoint of the ideal 
mT2 distribution in the squark-gluino production events is identical to the squark 
mass in our model points. 

A. 2 The effect of Matrix Element corrections to the signal distribution 

In this appendix we consider the matching effect of multi-jet matrix elements (ME) 
and parton showers on the distribution. When we calculate the signal 171x2 dis- 
tributions in this text, we generate SUSY processes at the lowest-order hard process 
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and then generate multi-jet events by parton showers. In general, there are ME 
corrections from hard parton emissions in the lowest order hard process, which may 
not be included in the parton shower approach. Note that we have applied cuts 
Pt > 50 GeV and |?7| < 3 for the jets to be included in the hemispheres. We need to 
check that this is enough to kill the effects of initial state radiations. 

When the ME corrections are taken into account, we should avoid double count- 
ing of emissions in overlapping phase space and need some kind of matching scheme 
to merge the ME corrections. Here, we study the ME corrections using the MadGraph 
/MadEvent MC generator [Q, in which the matching between the ME corrections 



and the parton showers is implemented. For our analysis, we use a modified MLM 
matching procedure with kj_ jets. In this scheme, the parton emissions are separated 
into two phase space regions at some kj_. In MadGraph/MadEvent, only events with 
enough separated partons, k± > xqcut, are generated after the matrix element sim- 
ulation. Then parton showering is performed and the partons are clustered into jets 
using the kj_ algorithm. After this procedure, the matching between the jets and the 
partons from the matrix elements is performed using Pythia. If the distance between 
them is larger than Qcut, the event is discarded in order to avoid double counting. 

In order to see the effect of the additional jet emission, we generate the SUSY 
events for the mSUGRA point SPS la using MadGraph/MadEvent. The generated par- 
ton level events are interfaced with Pythia to take into account the matching and the 
hadronization. We take the matching parameters as xqcut= 40 GeV, Qcut= 60 GeV. 
After hadronic events are generated, we use AcerDET for detector simulations. We 
apply the same cuts given in the Sec. 4 to select the events. 

In Fig. |l^(a), we plot the mT2 distributions for gluino pair-production processes 
with 0, 1, 2 jets. Since the total cross sections could receive large NLO corrections, 
the shape of distribution is more important. For comparison, we normalize the 
each distribution to unity. We can see that the shapes of the mT2 distributions are 
stable against the ME corrections. This is a good feature to obtain information on 
the gluino mass from the endpoint of the 771^2 distributions. 

In Fig. |l^(b) and 12(c), we plot the mT2 distributions for squark pair-production 
(squark-gluino) processes with 0, 1 jet. We also normalize each distribution to unity. 
Again we can see that the mT2 distributions are rather insensitive to the ME correc- 
tions and the matching. 
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